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Abstract 

D ■ 
£/} ■ In this work, we study six tree-dominated B — > ai(1260)ai(1260) and 6i(1235)&i(1235) decays 

in the perturbative QCD(pQCD) approach, where a\(bi) is a 3 Pi( 1 Pi) axial-vector meson. Based 

on the perturbative calculations and phenomenological analysis, we find that: (a) the CP- 

averaged branching ratio of B° — > a{a^ decay in the pQCD approach is 54.7 x 10 -6 , which 

agrees well with the current data and the predictions given in the QCD factorization approach 

Q-r within errors; (b) the numerical results for the decay rates of other five channels are found 

to be in the order of 10 -6 ~ 10 -5 , which could be accessed at B factories and Large Hadron 

Collider(LHC) experiments; (c) other physical observables such as polarization fractions and 

direct CP-violating asymmetries are also investigated with the pQCD approach in the present 

work and the predictions can be confronted with the relevant experiments in the near future; 

£f~) ■ (d) the different phenomenologies shown between B — > a\a\ and B — > bibi decays are expected 

to be tested by the ongoing LHC and forthcoming Super-B experiments, which could shed light 

on the typical QCD dynamics involved in these decay modes, as well as in 3 Pi meson a\ and 

1 P\ meson b\. 
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I. INTRODUCTION 

Charmless B meson decays to final states involving two axial- vector mesons (AA) have 
attracted attentions in theory and experiments in the last few years PHZ|. It is expected 
that through the study of B — > AA decays, the issues related with the internal structure, 
such as the angles between the mixtures of 3 Pi and/or 1 Pi states [8|, |9J], of the light 
axial-vector mesons can receive new understanding. On one hand, B — > AA will open 
another window to study their physical properties; The CP asymmetries of these decays, 
on the other hand, shall provide another way to measure the Cabibbo-Kobayashi-Maskawa 
(CKM) angles (3 and a. Furthermore, analogous to B — >• VV decays, being constructed 
of three polarization states, the charmless B decays to AA mesons are expected to have 
rich physics and provide much more information on the underlying helicity structure of 
the decay mechanism through polarization studies [3|. 

Very recently, the measurement on branching ratio(BR) and fraction of longitudinal 



polarization^) for B° — > a\ ( 1260) + a\( 1260) [5|, |l0) decay has been reported by BaBar 
Collaboration, 

BR^ -> a 1 (1260) + ai(1260)-) Exp . = 47.3 ± 12.2 x 10~ 6 ; (1) 

f L (B° ->• a 1 (1260) + ai(1260)") Ex p. = 0.31 ± 0.24 . (2) 

with large uncertainties. However, this measurement will be improved rapidly with the 
running of Large Hadron Collider(LHC) experiments, in which the events of B mesons are 
expected to be produced more than those collected in the B factories by about 3 orders 
per year. 

On the theory side, B — > aiajj decays have been studied in the literature [2|-l4|, but the 
predictions on BRs of the considered channels are significantly different from each othe r by 
employing the approach with naive factorization [ll| and QCD factorization(QCDF) |12| . 
respectively. For B° — > afa7 mode, for example, the branching ratio predicted in naive 
factorization is 6.4 x 10~ 6 [2j, while that presented in QCDF is 37.4 x 10~ 6 [3j]. One can 
easily see that the former result is too small to be confronted with the preliminary data 
and the latter one going beyond the naive factorization is large enough to be compatible 
with the measurements. As the counterparts of B — > aiOi decays, B — >■ b\b\ modes have 
also been investigated within the framework of QCDF [3] and the BRs are found to be 
in the order of 10~ 6 ~ 10~ 5 within large uncertainties. More important, an interesting 
pattern of the BRs for B — > b\b\ decays is exhibited through the calculations based on 
QCDF in Ref. Q, i.e., BR(B° -> b\b\) > BR(B + ->■ 6+6?) > BR(B° -> 6+6^), which is 
significantly contrary to that for B — > c^ai channels, i.e., BR(5° — > a^af) < BR(5 + — > 
afa®) < BR(£>° — > a^a^). Here, we want to mention that, as stated in Ref. [31], the 
troublesome endpoint singularities from hard spectator scattering and annihilation decay 
amplitudes always exist in the framework of QCDF and have to be determined through 
the input parameters fitted from the relevant precision measurements. 

Inspired by the above interesting facts from both theoretical and experimental aspects, 
we here study B — > aiOi and b\bi decays in the present work by employing the low energy 
effective Hamiltonian 13j and the perturbative QCD(pQCD) approach |14l. Il5l . 



1 Hereafter, for the sake of simplicity, we will adopt the forms a\ and 6i(to be shown below) to denote 
the meson ai(1260) and 6i(1235) in the content, respectively, unless otherwise stated. 



It is worth of stressing that the nonfactorizable spectator and annihilation diagrams 
are calculable perturbatively in the pQCD approach. Furthermore, the new measure- 



ments on the pure annihilation B s — > n + n decay reported by CDF 16 1 and LHCb [17 



collaborations last year confirmed the previous pQCD predictions |18H20J , while the mea- 
sured large decay rate of B° — > K + K~ also naturally explained by the renewed pQCD 
predictions [20| . 

Although <2i and 61 mesons embrace the same components at the quark level, because 
of different couplings of orbital and spin angular momenta, as explored in the QCD sum 
rule method [8], the hadron dynamics of &i is very different from that of its partner, 
a\. In our numerical evaluations, the different phenomenologies do appear between these 
considered B — y aidi and B — )■ b\bi decays. Therefore, one could expect reasonably 
that more new information on the contents such as polarizations, CP asymmetries, CKM 
unitary angles, even the knowledge of color-suppressed processes in B meson decays [2l|- 
24| may be deduced through the detailed studies on B — > a x ai and b±bi decays. Moreover, 
the hadron dynamics could be implicated by these perturbative calculations with the help 
of precision measurements. 

The paper is organized as follows. In Sec. [TTJ we present the theoretical framework 
on the low energy effective Hamiltonian, formalism of pQCD approach and mesons' wave 
functions. Then we perform the perturbative calculations for B — > a x ai and &1&1 decays 
in Sec. IIHI The analytic expressions of the decay amplitudes for the considered modes 
are also grouped in this section. The numerical results and phenomenological analysis are 
given in Sec. IIVI The main conclusions and a short summary are presented in the last 
section. 



II. THEORETICAL FRAMEWORK 

For the considered B — > qiaj and b±bi decays with b — > J transition, the related weak 



effective Hamiltonian H e g 13j can be written as 



H eS = 2| J v: b V ud [C x {^)O u M + CMO u M] - V: b V td JT Ci(/*)0 4 (/*) 1 , (3) 

with the Fermi constant Gf = 1.16639 x lCT 5 GeV~ 2 , CKM matrix elements V, and 
Wilson coefficients Ci(n) at the renormalization scale /i. The local four-quark operators 
Oi{i = 1, • • • , 10) are 

O" = (d a Ui3)v-A(u/3b a )v-A , Ol = (d a u a ) V - A (u/3bi3)v-A ; 

03 = (d a b a )v-A^2(q'pq'i3)v-A , 4 = {d a bp)v-A ^2((?/3<la)v-A , 

<?' g' 

5 = (dAlv-A^l^^+A! 6 = (dgb^y^A y ^q'pq'^y+A ; 



q' q' 

Q Q 

07 = -{dJj^y^A^eq^q'pq'^v+A , 8 = -{d a bp) V - A 2je g /(^^)y +A , 

q> q' 

Q Q 

09 = ■^{d ol b ol )v-A} j e q '{q' q' )y-A , 0io = 77 {d a bp)v- A }, e qf (^ q' a )v- A 

q> q> 



(4) 



with the color indices a, (3 and the notations (q'q')v±A = 9'Tm(1 ± 7s)?'- The index q' in 
the summation of the above operators runs through u, d, s, c, and b. Since we work in 
the leading order[C(a s )] of the pQCD approach, it is consistent to use the leading order 
Wilson coefficients. For the renormalization group evolution of the Wilson coefficients 
from higher scale to lower scale, we use the formulas as given in Ref. [14| directly. 

The pQCD approach, one of the method based on QCD dynamics in the market, has 
been employed to treat two-body nonleptonic B^ decays extensively. As a unique feature 
different from other two factorization approaches, i.e., QCDF and SCET (soft-collinear 
effective theory) [25], the pQCD approach is based on the framework of kx factorization 
theorem with taking the tranverse momentum k?, generally considered as a small and 
negligible scale, of the valence quarks in the hadrons into account, which results in the 
Sudakov factor smearing the endpoint singularities in the decay amplitude and makes the 
nonfactorizable spectator and annihilation diagrams perturbatively calculable, aside from 
the emission one. 

Because of the rather heavy b quark, for convenience, we will work in the rest frame 
of B meson. Throughout this paper, we will use light-cone coordinate (P + ,P~,P T ) to 
describe the meson's momenta with the definitions 

p± = Po±Ps and p T = ( piiP2 ). (5) 

Then for B° — > a±a± decay, for example, the involved three meson momenta in the 
light-cone coordinates can be written as 

P 1 = ^(1,1,0 T ), P 2 = ^(l-r 3 2 ,r 2 2 ,0 T ), P 3 = ™|( r 2,l- r 2,0 T ), (6) 

respectively, where the af (aj~ ) meson moves in the plus (minus) z direction carrying the 
momentum P 2 (P3) and r 2 = r 3 = m ax jvfiB- The longitudinal and transverse polarization 
vectors of axial-vector meson are denoted by e L and e T , respectively, satisfying P ■ e = 
in each polarization. The longitudinal polarization vectors, e 2 and ef , can be chosen as 

* mB (l-r 3 2 ,-r 2 2 ,0 T ) and e L 3 = -^(-r 2 , 1 - r 2 , T ). (7) 



V2m ai " ~ " \pl 



m n 



And the transverse ones are parameterized as e^ = (0, 0, It) and ej = (0, 0, It). 

Putting the (light) quark momenta in B, af and a{ mesons as fci, k 2 , and k 3 , respec- 
tively, we define 

k 1 = (x 1 P+,0,k 1T ), A; 2 = (x 2 P 2 + ,0,k 2T ), fc3 = (0,X3P 3 -,k3T). (8) 

Then, for B° — > a^a^ decay, the integration over fcf , k%, and k£ will conceptually lead 
to the decay amplitude in the pQCD approach, 



A(B -)• afa 1 ) ~ /( 



dxidx2dx'ibidbib2db2b'idb'i 



•Tr 



C(t)^ B (x 1 ,b l )^ at (x 2 ,b 2 )^ a -(x 3 ,b 3 )H(x h b i ,t)S t (x l )e-^ . (9) 



where b{ is the conjugate space coordinate of k iT , and t is the largest energy scale in func- 
tion H(xi,bi,t). The large logarithms ln(mw/t) are included in the Wilson coefficients 



C(t). The large double logarithms (In x*) given rise from loop corrections to the weak 



decay vertex are summed by the threshold resummation 26j, and they lead to S t {Xi) 
which can decrease faster than any power of x as x — > 0, then remove the endpoint 
singularities. The last term, e~ s ^\ is the Sudakov factor which suppresses the soft dy- 
namics effectively [27|. Thus it makes the perturbative calculation of the hard part H 
applicable at intermediate scale, i.e., m^ scale. We will calculate analytically the func- 
tion H(xi,bi,t) for the considered decays at leading order in a s expansion and give the 
convoluted amplitudes in next section. 

The pQCD predictions depend on the inputs for the nonperturbative parameters such 
as the decay constants and distribution amplitudes. For heavy B meson, in principle, 
both Lorentz structures of the wave function should be considered in the calculations. 
However, the contribution induced by the second Lorentz structure is numerically small 
and approximately negligible 28] , we therefore employ the following set of heavy B meson 
wave function [l4J . 



<$> B (x,b) 



y/E 



) +m B )j 5 (p B (x,b) 



where the distribution amplitude 4> B (x,b) has been modeled as [141 ]. 



4>b{x, b) = N B x 2 (l — x) 2 exp 



2 1,2 



1 ( xm B \ u z B b 



Ufi / 



(10) 



;n) 



In recent years, the shape parameter u B in Eq. ffTTj) has been fixed at 0.40 GeV in the 
pQCD approach by using the rich experimental data on the B mesons with f B = 0.19 GeV. 
The normalization factor N B is related to the decay constant f B through 



/ 



B 



dx(j) B (x,b = 0) 
o 2a/d 



(12) 



Correspondingly, the normalization constant N B is 91.745 for u B = 0.40. To analyze 
the uncertainties of theoretical predictions induced by the inputs, we will vary the shape 
parameter u B by 10%. 

For the wave functions of axial-vector oi and b\ mesons, one longitudinal (L) and two 
transverse(T) polarizations are involved, and can be written as [29j , 

1 

y/E 



*a(x) 



-^lA m A i f<f> A {x) + i A L f<p A (x) + m A <p s A (x) 



a0 



$ T A (x) = -f=iAm A { 



■ *T iv 



A " 'A\ X ) + F A 



*T) 



b T A (x) + m A ie ilvpal ^^n p v u <\) a A {x) 



(13) 
(14) 



where x denotes the momentum fraction carried by quark in the meson, and n = (1, 0, Oy) 
and v = (0, 1, 0<r) are dimensionless light-like unit vectors. We here adopt the convention 
e 0123 = 1 for the Levi-Civita tensor e^^ '. 

The twist-2 distribution amplitudes for the longitudinally and trasversely polarized 
axial- vector 3 P 1 and 1 Pi mesons can be parameterized as [a. l29| : 

3/a 



cj) A (x) 



J A 



X 



3/4 

V2Nc 



x(l — x) 
x(l — x) 



l 0A + 3fl lA ( 2x ~ !) + a 2A 2^ 2X ~ X ) Z 



3 
3 



1) 



i 0A + 3a 1A (2x - 1) + a 2A -(5(2x 



(15) 
(16) 



where f A is the decay constant. Here, the definition of these distribution amplitudes 
(J>a(x) an d <J>a(. x ) satisfy the following relations: 



(pa p. (X) = ; 

1 



iT I \ _ _!_ J Z P\ 

PzPiK 00 ) — a o 3 Pi; 



2^/2N c 



n Pl {x) = a 



h 



Pi 



1 Pi 



2^/2N c 



:/' , .s _ / J A 



0i Pl (x) 



2v^iVc 



(17) 



where a{J 3p = 1 and ajjip = 1 have been used. 

As for twist-3 distribution amplitudes for axial-vector meson, we use the following 
forms 29 : 



(f> f A (x) 
(f> s A (x) 



^= Li A {2x - l) 2 + \ a\ A (2x - l)(3(2x - l) 2 - 1) 



3/a d 

2y/2~W c dx 



{x(l-x)(da A + ai A (2x-l))} ; 



(18) 
(19) 



<j>\{x) 
d> a A (x) 



3// 



4V2N, 

Va d 
4y/2N~ c dx 



1 || 
2 %A 



1 + (2x - l) 2 ) + a\ A (2x - l) 3 
x(l - x)(a|) J 4 + a\ A (2x - 1)) 



(20) 
(21) 



The Gegenbauer moments a- A have been studied extensively in the literatures (see 
Ref. |8j and references therein), here we adopt the following values: 



f ai =0.238 ±0.010 GeV, 
f bl =0.180 ±0.008 GeV, 



2,ai 



'!,«* 



-0.02 ±0.02 
-1.95 ±0.35 



"l,oi 



l 2,6i 



-1.04 ±0.34; 
0.03 ±0.19. (22) 



Note that we have included the intrinsic b dependence for the heavy meson wave 
function <pB but not for the light axial- vector meson wave function <p A . It has been shown 
that the intrinsic b dependence of the light meson wave functions is not important and 
negligible 30j, |3l|. It is reasonable to preliminarily assume that the intrinsic b dependence 
of the Oi and b\ wave functions, which are still unknown, is not essential either. 



III. PERTURBATIVE CALCULATIONS IN PQCD APPROACH 



There are three kinds of polarizations of a axial-vector meson, namely, lon gitu dinal 
(L), normal (N), and transverse (T). Analogous to the B — >■ pp decays J32I434J . the 
amplitudes for the B — >■ aiai decays are also characterized by the polarization states of 
these axial- vector mesons. In terms of helicities, the decay amplitudes Ai^ for B — > 
ai(P 2 , e2) a i(-^3' e 3) decays can be generally described by 



M (a) 



e 2li [cr)€i,.[cr) 



-3v\ 



ag*" + 



b c 

-P?P? + i- 



s^P 2a P w 



m 2 B M L + m 2 B M N e* 2 (a = T) ■ e* 3 (a = T) 
+iM T e a ^e* 2a (a)e$p(a)P 27 P 3p , 



(23) 



( 



-f -f 

V V 

£"•"* 1 i s ~~8 1 A" 

I I"; e | I I", e | 



fa; 



ft) 





(9) (h) 



FIG. 1. (Color online) Typical Feynman diagrams for B — > a\a\ decays at the lowest order in 
the pQCD approach. By replacing the 3 Pi meson a± in (a)-(h) with the 1 P\ meson b\, one will 
obtain the corresponding Feynman diagrams for B — >• b\b\ decay modes. 

where the superscript a denotes the helicity states of two mesons with L(T) standing for 
the longitudinal (transverse) component. And the definitions of the amplitudes Aii(i = 
L, N, T) in terms of the Lorentz-invariant amplitudes a, b and c are 



m 2 B M L = a e*(L) • e^(L) + 



m ai m ai 



-e* 2 (L) ■ P 3 e*(L) • P 2 



m% M 



N = a 



(24) 



m\ M T 



r 2 r 3 



We therefore will evaluate the helicity amplitudes M.L-, -M-N, -M.T based on the pQCD 
approach, respectively. 

As illustrated in Fig. [Q there are 8 types of diagrams contributing to the B —¥ a±ai 
and &1&1 decays at the lowest order in the pQCD approach. We firstly calculate the usual 
factorizable spectator (fs) diagrams (a) and (b), in which one can factor out the form 
factors B — > a\ and B — >■ b\. The corresponding Feynman amplitudes with longitudinal 
polarization(L) are given as follows, 

(i) (V - A)(V - A) operators: 



F L 



87rCjpf ai m% / dxidxz / bidbib 3 db 3 
Jo Jo 

x0 B (x!, 60 { [(1 + x 3 )(j) ai (x 3 ) + r ai (1 - 2x 3 )« (x 3 ) + < (a*))] 
xh fs (x 1 ,x 3 ,bi,b 3 )E fs (t a ) + 2r ai 4> s ai (x 3 ) h fs (x 3 ,x 1 ,b 3 M) E fs (t b )} , (25) 



where Cp = 4/3 is a color factor. The convolution functions Ej, the factorization 



hard scales ti, and the hard functions hi can be referred to Ref. [35 



(ii) {V - A)(V + A) operators: 



Fg P1 = -Fk , (26) 



fs - ± fs ' 

which is originated from (ai|V + A\0) = — (ai\V — A\0). 
(iii) (S — P)(S + P) operators: 

Fff 2 = ; (27) 

because the emitted axial-vector meson can not be produced through a scalar or a 
pseudoscalar current. 

For the nonfactorizable spectator (n/s) diagrams 1(c) and 1(d), the corresponding de- 
cay amplitudes can be read as 



(i) (V - A)(V - A) operators: 

32 f 1 f°° 

M^f s = —=ixCFm 2 B I dx\dx 2 dx 3 / bidbib 2 db 2 4>B(xi,bi) 

V6 Jo Jo 

X(f) ai (x 2 ) {[(1 - x 2 )(f) ai (x 3 ) - r ai x 3 ((f) ai (x 3 ) -^(xa))] 

xE nfs (t c )h c nfs (x 1 ,x 2 , x 3 , 6i, b 2 ) - h* fa (xi,X2, x 3 , h, b 2 ) 

x [(x 2 + x 3 )(f> ai (x 3 ) - r ai x 3 ((f) s ai (x 3 ) + ^(xs))] E nfs (t d )} , (28) 

(ii) (V -A)(V + A) operators: 

32 r 1 C 00 

M nU = —=^Cpm 1 B \ dx\dx 2 dx 3 \ bidb\b 2 db 2 (f)B{x\,bi) 
V6 Jo Jo 

xr ai { [(1 - x 2 )((f) s ai (x 2 ) + 0* ai (x 2 ))0 ai (x 3 ) - r ai {<f> s ai (x 2 ) 
x [(x 2 -x 3 - l)(p s ai (x 3 ) - (x 2 + x 3 - 1)0^ (x 3 )} + <f) ai (x 2 ) 
x[(x 2 + x 3 - l)(j) s ai (x 3 ) + (1 - x 2 + x-i)^])] h c nfs (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 

xE nfs (t c ) - h^ lfs (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )E nfs (t d ) [x 2 ((f) s ai (x 2 ) - (f) ai {x 2 )) 
X(f) ai (x 3 ) +r ai (x 2 ((j) ai (x 2 ) - 4> ai (x 2 )) ((f) ai (x 3 ) - <0r 3 )) 

+x 3 (<p s ai (x 2 ) + <(x 2 ))«(x 3 ) + <(x 3 )))]} , (29) 

(iii) (S - P)(iS + P) operators: 



-*1 /»oo 

M^f s z - —=nCFm B I dx\dx 2 dx 3 I bidbib 2 db 2 (f>B(xi,bi] 



f L;P2 32 _/ -, ...2 



o ./o 



X0 ai (^2) {[(^2-^3- l)0oi(aJ3) + r ai X 3 {(t) ai {x 3 ) +<t> ai {x 3 ))] 

xE nfs (t c )h c nfs (xi, x 2 , x 3 , bi, b 2 ) + h^ lfs (x 1 ,x 2 , x 3 ,h,b 2 ) 

x [x 2 <p ai (x 3 ) - r ai x 3 ((f) ai (x 3 ) - 4> ai (x 3 ))] E nfs (t d ) } ; (30) 



In the above three formulas, i.e., Eqs. ( |28l) -( l30j) . one can find that there exist cancellations 
between the contributions of the two diagrams in Fig. [T](c) and[TJ(d). 

The Feynman diagrams shown in Fig.[Tfe) and[Tff) are the nonfactorizable annihilation(n/a) 
ones, whose contributions are 

(i) (V -A)(V- A) operators: 

'1 />oo 



>2 



32 f f°° 

M nfa = —^iiCpm 2 ^ \ dx\dx 2 dx 3 / bidbib 2 db 2 
V6 Jo Jo 

X(p B (x 1 ,b 1 ) {[(I - x 3 )(j) ai (x 2 )(j)a 1 (x3) + r ai r ai (^M 
x [(1 + x 2 - x 3 )(f) s ai (x 3 ) - (1 - x 2 - z 3 )0* ai ( x s)] + < (x 2 ) 
x [(1 -x 2 - x 3 )<p s ai {x 3 ) - (1 + x 2 - £3)0*^X3)])] £ n / a (*e) 
x K fa ( Xl , x 2 , x 3 , h, b 2 ) - E nfa {t f )h{ fa {^^ xs, b u b 2 ) 
x [x 2 (j) ai (x 2 )<f>ai M + r ai r ai {<p s ai (x 2 ) [(x 2 -x 3 + 3)0^ (x 3 ) 
+(1 -x 2 - x 3 )0* ai (x 3 )] + 0* 1 (a ; 2)[(^2 + x 3 - l^fo) 
+(l-x 2 + x 3 )<(x 3 )])]} , (31) 

(ii) (V -A)(V+ A) operators: 

M ' f = —j='KCFvr? B \ dx\dx 2 dx 3 / b\db^b 2 db 2 
V6 Jo Jo 

x0 B (xi,6i) {[r ai x 2 (C 1 (x 2 ) + 0* ai (x 2 )) fll (x 3 ) -r ai (l - x 3 ) 

X 00! (^2) (0a! (^3) -^ai^s))] E n fa(t e )K fa (x l ,X 2 ,X 3 ,b l M) 

+ [r ai (2-x 2 )(0 ai (x 2 ) + ai (x 2 ))0 ai (x 3 ) -r ai (l + x 3 ) 
x0a 1 (a; 2 )(0: i (a; 3 )-0 ai (x 3 ))] J E n/a (t / )^ /a (x 1 ,x 2 ,x 3 ,6 1 ,6 2 )} , (32) 

(iii) (S — P)(S + -P) operators: 

M n ' fa 2 = -=7rC F m 2 B / dx\dx 2 dx 3 \ b\db\b 2 db 2 

V6 Jo Jo 

X(j) B (xi,bi) {[(1 -x 3 )0 ai (a; 2 )0 ai (a; 3 ) +r ai r ai (0 ai (x 2 ) 
x [(x 2 - x 3 + 3)0^ (x 3 ) - (1 - x 2 - x 3 )0 ai (x 3 )] + ai (x 2 ) 
x [(1 - x 2 - x 3 )0 ai (x 3 ) + (1 - x 2 + x 3 )0* ol (x 3 )])] E nfa (t f ) 
x Kfa( x i> x z> x 3, 61, h) - E nSa {t^h e nfa {x u x 2 , X 3 , 61, 62) 
x [x 2 ai (x 2 )0 ai (x 3 ) + r ai r ai (0^ (x 2 ) [(1 + x 2 - x 3 )(j) s ai (x 3 ) 

+ (1 - X 2 - X 3 )0 ai (x 3 )] + Q1 (X 2 )[(X 2 + X 3 - l)0 ai (x 3 ) 

-(l + x 2 -x 3 )0 ai (x 3 )])]} ; (33) 

For the last two diagrams in Fig. [TJ i.e., the factorizable annihilation(/a) diagrams 
[TJg) and[T](h), we have 

9 



(i) {V - A)(V - A) operators: 



/*1 /*oo 

Ff a = -8irC F m 2 B / dx 2 dx 3 / b 2 db 2 b 3 db 3 {[x 2 (p ai (x 2 )(f) ai (x 3 ) + 2r ai r ai 
Jo Jo 

X ((X2 + 1)C (X 2 ) + (X2 - 1)< (Z2)) C (X 3 )] M^, 1 - X 3 , 6 2 , 63) 

xE fa (tg) - [(1 - x 3 )0 ai (x 2 )0 ai (a; 3 ) - 2r ai r ai 0' 1 (x 2 ) ((z 3 - 2)0*, (x 3 ) 

-X 3 0^(X 3 ))] ^/ a (tfc)/i/ a (l-X 3 , ^2,63,62)} , (34) 



(ii) (V - A)(V + A) operators: 

'fa - t'fa ' 



*£" = Fk » (35) 



(iii) (5 — P)(*S + P) operators: 

F^'f 2 = lQ-nC F m 2 B j dx 2 dx 3 / b 2 db 2 b 3 db 3 {[2r ai (l>a 1 (x 2 )(f) s ai (^a) 
Jo Jo 

+r oi x 2 (0a 1 (a;2) - ^(a^))^^)] h fa {x 2 ,l - x 3 ,b 2 ,b 3 )E fa (t g ) 
+ [2r ai C x (x 2 )0 ai (rr 3 ) + r ai (1 - x 3 )<f) ai {x 2 ){4> s ai (x 3 ) + < (x 3 ))] 
xEf a (t h )hfa(l - x 3 ,x 2 ,b 3 ,b 2 )} . (36) 

It is interesting to notice that there is a large cancellation in the Fk, i.e., Eq. (I34j) . from 
the factorizable annihilation diagrams Hfg) and H^h) , which can result in the exact zero 
contribution in the SU(3) limit. 

We can also present the factorization formulas for the Feynman amplitudes with trans- 
verse polarizations, 

Fj? s = 87rC F f ai m 2 B / dx x dx 3 \ bxdbxb 3 db 3 (f) B (xx,bx)r ai 
Jo Jo 

X { K,(x 3 ) + r ai x 3 (KM - 0^(x 3 )) + 2r ai # 1 (x 3 )] E fs (t a ) 

x/i /s (xi,x 3 ,6i,& 3 ) +r ai (0^(x 3 ) +<f>l 1 (x 3 )) h fs (x 3 ,Xx,b 3 ,bi) E fs {t b )} , (37) 

Fj s = l£mC F f ai m 2 B \ dx x dx 3 \ b 1 db 1 b 3 db s <f> B (xi : bx)ra 1 

Jo Jo 

x { K(z 3 ) - r ai x 3 {<t> v ai {x 3 ) - r ai (xs)) + *r ai (f> a ai M] E fs (t a ) 
xh fs (xx,x 3 ,bx,b 3 ) + r ai ((j) a ai (x 3 ) + (j) v ai (x 3 )) h fs (x 3 ,Xx, h, bx) E fs {t h )} , (38) 

(39) 
(40) 

(41) 
(42) 



p JV;Pi 


J7 N 


b u ~ 


-Fjs ; 


b fs 


= 0, 


T?T;P 2 
b fs 


= 0; 
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32 f 1 f°° 

M^r s = —=TcC F m 2 B / dx x dx 2 dx 3 / bidbib 2 db 2 4>B(xi,bi)r ai 

v6 Jo Jo 

x { [(1 - x 2 )(<j) a ai (x 2 ) + <f) v ai (x 2 ))^i fas)] ^/ s (^i, ^2, a*, W, 62) 

xE nfs (t c ) +hi f9 (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) [x 2 (<p a ai (x 2 ) + <f ai (x 2 ))4> T ai (x 3 ) 

-2r ai (x 2 + x 3 )(r ai (x 2 )r ai (x 3 ) + <j? ai (x^l (x s ))] £„/. ft,) } , (43) 

f\A r^~ /*oo 

ML = —=ixC F m 2 B I dxidx 2 dx 3 / bidbib 2 db 2 (f> B (xi,bi)r ai 
V6 Jo Jo 

X { [(1 - X 2 )(0ai (^2) + € x (x2)W ai M] Kfsi^l, X 2 , X 3 , h, 6 2 ) 
Xf? n/a (t c )+/lj| /g (xi,X2,X3,6l,62) [^(^ (X 2 ) +0a 1 (^2))0l 1 M 

-2r ai (x 2 + a*) « W< (*s) + < (^)C (*s))] £„/. (t„) } , (44) 



32 Z" 1 Z" 00 

Mnfs 1 =: 7^C F rr? B \ dx 1 dx 2 dx 3 / b 1 db 1 b 2 db 2 (f) B (x 1 ,b 1 )r ai 

Jo Jo 



nfs ^ 



*X 3 4> T ai {x 2 ){<p a ai {x 3 ) -4> V ai {x 3 )) [^ /a (xi,X 2 ,X3,6i,6 2 )-Bn/ a (*c) 

+/i^/ a (a:i,a;2,a:3,6i,62)^n/»(td)] , (45) 

M™=2M?>?; (46) 

M nu 2 = i=nC F m 2 B dx 1 dx 2 dx 3 / b 1 db 1 b 2 db 2 (f) B (x 1 ,b 1 )r ai 

V6 Jo Jo 

x {[^2(0^(^2) -^oi^a))^ (a*)] /ij| /s (xi,x 2 , 0:3,61,62) 
xE nfs (t d ) + h c nfs (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) [x 2 (<j) a ai (x 2 ) - <t> v ai {x2)W ai {x 3 ) 
+2r ai (l - x 2 + ^(^(^^(xa) - ft^Cfe))] £»/.&)} , (47) 

/^/| /*1 /*oo 

M^f 2 = —irC F m B dxidx 2 dx 3 \ b 1 dbib 2 db 2 (f) B (xi,b 1 )r ai 

V6 Jo Jo 

X (M^O^) -<(l2))ft(l3)] ^ /g (xi,X 2 , ^3,61,62) 

xE nfs (t d ) +h c nfs (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) [x 2 {(p a ai {x 2 ) - (i) v ai {x 2 ))(t) T ai (x 3 ) 
+2r ai (l - x 2 + X3)(C 1 (x 2 )^ 1 (x 3 ) - CO^Cfo))] ^n/ s (t c )} , (48) 

f\A r^ /*oo 

M nfa = i=nC F m 2 B dx 1 dx 2 dx 3 / bidbib 2 db 2 r ai r ai 

v6 Jo Jo 

x [«> 2 )Cte) + <(*W] I • ( 49 > 

1 28 Z" 1 /"°° 

M nfa = i=nC F m 2 B / dx 1 dx 2 dx 3 / b x db x b 2 db 2 r ai r ai 

V6 Jo Jo 

x [C, foX (*s) + < (* 2 )C (a*)] } , (50) 
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\/6 ,/o ./o 



M„c'o = —^TrOpinr, / dx\dx 2 dx 3 / b\db\b 2 db 2 

} VQ Jo Jo " " ~ " 

x0 B (xiA) { [^1^2(00! (^2) + 0a 1 (^2))0T 1 (s 3 ) -r ai (l -s 3 ) 



x<(^2)«(s 3 ) -0^(s 3 ))] E nfa (t e )h e nfa (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 
+ [r ai (2 - s 2 )(0^(s 2 ) + C 1 (^))0r i (s 3 ) - r ai (l + s 3 ) 
x0r i (^2)(C 1 (s 3 )-C 1 (s 3 ))]K/a(i/)^ /a (si,s 2 ,s 3 ,6 1 ,6 2 )} , (51) 



M^=2M^; (52) 

<f = <a , (53) 

^5r 2 = -Kfa ; (54) 

Ff a = 87rC F m 2 B / dx 2 rfs 3 / b 2 db 2 b 3 db 3 r ai r ai { [4> a ai {x 2 ) ((s 2 + 1)0^ (s 3 ) 
Jo Jo 

+(s 2 - 1)^(^3)) + C(s 2 ) ((s 2 + l)0" ai (a*) + (s 2 - 1)0^ (s 3 ))] £/.(t„) 

X/l /a (s 2 , 1 - S 3 , 62, 63) + [(*3 - 2)(0^(S 2 )C(S 3 ) + 0^(S 2 )0^(S 3 )) 
-S 3 (C 1 (S 2 )0^ 1 (S 3 ) +^ 1 (X 2 )^ 1 (X 3 ))] ^/a(*h)/l/o(l-X3,X 2 ,6 3 ,62)} , (55) 

pi /'OO 

Ff a = 16irC F m B dx 2 dx 3 / b 2 db 2 b 3 db 3 r ai r ai { [4> v ai {x 2 ) ((s 2 + 1)0^ (s 3 ) 
Jo Jo 

+ (s 2 - 1)0^ (s 3 )) + 0^(s 2 ) ((s 2 + l)0 u ai (s 3 ) + (s 2 - l)0^(s 3 ))] E fa (t a ) 
x fc /a (a; 2 , 1 - s 3 , 6 2 , 6 3 ) + [(s 3 - 2) « (s 2 )0^ (s 3 ) + 0^ (s 2 )0^ (s 3 )) 
-s 3 (0^(s 2 )0^(s 3 ) + 0^(s 2 )0: i (s 3 ))] E fa (t h )h fa (l -x 3 ,x 2 ,b 3 ,b 2 )} , (56) 



*#* = *£ , (57) 

*£* = *£ , (58) 



pi poo 

F^ P2 = 16nC F m 2 B / dx 2 dx 3 / b 2 db 2 b 3 db 3 r ai 
Jo Jo 

X {K(^2) (0a 1 (^3) -0a 1 (s 3 ))] E fa (t g )h fa (x 2: 1 - S3, 62, 63) 

+ [«(^ 2 ) + 0: i (s 2 ))0j(x 3 ))] E fa (t h )h fa (l -s 3 ,s 2 ,fe 3 ,6 2 )} , (59) 

Fjf 2 =2F% p >. (60) 

Thus, for these considered six tree-dominated decay channels, by combining all the 
possible contributions from different Feynman diagrams, we can display the physical decay 

12 



amplitudes with three polarizations h = L,N,T as follows, 



M h (B° -> aUi) = A u 



A, 



(o 4 + Oio)^ 
1 



ai F% + CX /8 + C 2 M n h /a + oa/jjift 

+(C 3 + G,)M„% + (C 5 + C 7 )M n ^ + (C 3 + 2C 4 - - (C 9 - C 10 )) 
xM n ft /a + (C B - iCr)^ + (2C 6 + \C % )M%Z + (2a 3 + a 4 
+2a 5 + ha 7 + a 9 -a 10 ))f B F} a + (a 6 -^a s )f B F^ , (61) 



V2M h (B + ->■ a+a? 



) = A t 



( ai + a 2 )F f h s + (d + C 2 )M, 



n/s 



A, 



(2(C 9 + C; 



10 J 



-~(3C 7 + C 8 ))F£ + ~(C 9 + C w )M£ fa + *-C 7 M h n £ 



+ "C 8 M^ 



(62) 



V2M h (B° -> oSo?) = A u 



aa(/ B ^ - Ft) + d(M n V - M* 



-A t 



(04 - -(3a 7 



+3a 9 + 010))^, + (C 3 - ~(C 9 + 3Ci ))M^ - (C 5 - ^C 7 ) 



xM 



fe;Pi 



7-^2 



1 



- ~C S M^ + (C 3 + 2C 4 - -(C B - C 10 ))M 5 



2 ^)M n ^ 1 + (2C 6 + ^ 



ft 

n/a 



+ (C 5 - \c 7 )M h n f: + (2C 6 + \c s )M h n f: + (2a 3 + a 4 + 2a 5 



+^(a? - a 9 + a 10 ))f B F} a + (a 6 - -a s )f B F^ 



(63) 



In the above Eqs. (I6T]) - (I6"3"]) . A u and A* stand for the products of CKM matrix elements 
VubVud and V t * b V t d, respectively. The standard combinations Oj(r.h.s.) of Wilson coefficients 
are defined as follows, 



n 1 ^1 

ai = C 2 + — 



a-2 = Ci + — , 



Ci + —(i = 3- 10) . (64) 



where the upper (lower) sign applies, when i is odd(even). While for B —¥ bib± decay 
channels, one can easily obtain the analytic formulas for various decay amplitudes just by 
replacing ai with b\ in Eqs. (125]) - (163]) correspondingly. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will calculate numerically the CP-averaged BRs, polarization frac- 
tions, direct CP- violating asymmetries, and relative phases for those considered B — > oiOi 
and b\bi decay modes. In numerical calculations, central values of the input parameters 
will be used implicitly unless otherwise stated. 
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TABLE I. The CP-averaged predictions of physical observables for B° — > aja±,b~[bi decays 
obtained in the pQCD approach(This work). For comparison, we also cite the available experi- 
mental measurements [5j] and the theoretical estimates in the framework of QCD factorization [3|]. 





Decay Channels 


S° -> g+g, 


B° -> bfb. 


Parameter 


Definition 


This work 


QCDF 


Experiment 


This work 


QCDF 


Experiment 


BR(10~ 6 ) 


r/r tota i 


,. 7 +19.4+29.4+5.7 
OT - '-16. 9-23. 2-4. 4 


374 +16.1 + 9.7 
°'-— 13.7-1.4 


, 7 ,+ 10.5 + 6.3 
^' ■ -10.5-6.3 


,, ,+5.7+18.1 + 2.1 
^- -5.3-11.3-1.4 


, ,,+ 1.6+15.7 
11 -0.7-0.3 


- 


















fL 

f\\ 
f± 


\-A L \ 2 

I-AjJ 2 


, TK +0. 01+0. 03+0. 00 
u -' D -0. 00-0. 04-0. 00 
n , ,+0.00+0.02 + 0.00 
u,1 *-0. 00-0. 02-0. 00 
n in +0. 00+0. 02 + 0. 00 
u,1 "-0. 00-0. 01-0. 00 


64+ - 07 
U ' D *-0.17 


„ ,.+0.22 + 0.10 
u.o±_g 22-0 .10 


n qo + 0. 02 + 0. 04+0. 01 
u, °°-0. 01-0. 05-0. 00 
„ n7 +0. 01 + 0. 03 + 0. 00 
u,u ' -0.01-0.02-0.00 

n n ,-+o. oi+o. 02+0. oo 

u,u -0.01-0.01-0.00 


96+ 003 


- 


















0||(rad) 

0^(rad) 

A0||(rad) 
A^>j_(rad) 


n + arg -jJ-L 
T + arg -jj£ 


,, 1K +0. 01 + 0. 06 + 0. 01 
o,1D -0. 01-0. 04 -0.01 
„ 17 +0. 01 + 0. 06 + 0. 03 
°'-0. 00-0. 04-0. 01 

n n .+0. 00+0. 01+0. 00 

u,u -0. 01-0. 01-0. 01 
. R ,+0. 01 + 0. 01 + 0. 00 
i,D -0. 01-0. 01-0. 01 


" 


" 


9 r ,+0. 07+0. 07+0. 01 
z ' 01 -0. 06-0. 06-0. 01 
9 ,-+0.06+0.08 + 0.01 
z '*'-0. 05-0. 08-0. 01 

n .0+0.02+0. 07+0. 02 
u -^_ 0.04- 0.09- 0.04 
. ., + 0.04+0.09 + 0.02 
1 - lo -0. 02-0. 09-0. 03 


" 


- 


2 
0_L — 0J_ — ^ 




2 




■^CP 

.dir.L 
•^CP 

^CP 

,dir,J_ 

■"■cp 


r-r 
_r+r 

/ ir / i 


n ,,,+0.01 + 0.01 + 0.00 
u - u *_ 0.01 -0.01 -0.00 
n ,,+0.02 + 0.01 + 0.00 
u,± -0.02-0. 01-0. 01 

n (-9+0.08+0. 10+0. 04 

u ' JZ -0. 07-0. 07-0. 02 

n r ,+0.08 + 0.11 + 0.04 
u,o -0. 06-0. 06-0. 01 


" 


" 


n ,1,1+0. oi+o. oo+o. oo 

u - uu — 0.04— 0.01 — 0.00 

n ik+O.04+0.02+0.00 
uuo - 0.04- 0.03- 0.00 

n ,,+0.05 + 0.02 + 0.01 
u -^_ 0.06- 0.04- 0.03 

n c.o + 0.06 + 0.03 + 0.02 
U -^ B - 0.06- 0.04- 0.02 


" 


- 


/|i+/n 

f±-f± 
/, +/, 


- 



The QCD scale (GeV), masses (GeV), and B meson lifetime(ps) are 14i 136 



A 



MS 

m B 



0.250 

5.279 



m w 
m b 



80.41 

4.8, 



m 



a\ 



T B + 



1.23, 
1.638 



m bl 



1.21; 



t b o = 1.53 . 



(65) 



For the CKM matrix elements, we adopt the Wolfenstein parametrization and the 
updated parameters A = 0.814, A = 0.2257, p = 0.135, and fj = 0.349 . 

Utilizing the above chosen distribution amplitudes and the central values of the relevant 
input parameters, the resultant B — Y oi and B — > b\ form factors at maximal recoil, 



V f 







Hl ± 



V, 



B^b! 







34+ 010 

u -^ d -0.09 ' 



A B^ bl 



n 97+O.O6 
u ' z '-0.05 ' 



0.20 



+0.05 
0.04 5 



Vi 



B—>ai 



B^b x 



41+ - 10 ■ 

u -^ ± -0.08 ) 

™ +0 - 07 

u - ou -0.06 • 



(66) 
(67) 



associated with the longitudinal, parallel, and perpendicular components of the B — > q-ict i 
and B — > b\bi decays, respectively, are in good consistency with those as given in Ref. [29J. 
As a comparison, we quote the form factors used in the B — > aidi and b\bi decays in 
the QCD factorization [3j, 



t/B->«i 
v 



0.30 ±0.05 
0.39 ±0.07 






0.30 ±0.05 
0.16 ±0.03 



Vi 



B—^ai 



B^6i 



= 0.60 ±0.11; (68) 
0.32 ±0.06. (69) 



One can find that the form factors in the light-cone sum rule(LCSR), Eq. f)68|69l) . are 
basically consistent with those in the pQCD approach, Eq. (I66|67p . within errors. But, 
it should be noticed that, for B — >■ Oi transition, Vo ; lcsr < Vo ;p qcd in the longitudinal 
polarization while A L csr > ^4pQCD an d V^lcsr > K ; pQCD m both tranverse polariza- 
tions, which may reduce the polarization fraction fi, for example, for B° — > a^a^ mode 
evidently. 
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A. CP-averaged Branching Ratios 



For B — y O1O1 and b\ &i decays, the decay rate can be written explicitly as, 

Y^ M {a)] M {u) 



r 



167nn 2 B 



(70) 



=L,T 



where |P C | = |P2z| — |P3z| is the momentum of either of the outgoing axial- vector mesons 
and M {a) can be found in Eqs. (|6B-(|631). 

The theoretical predictions on CP-averaged BRs for B — > aidi and b\ &i decays evalu- 
ated in the pQCD approach, together with the results in QCDF approach and the available 
experimental data, have been grouped in Tables Hl lIIII The first error is from the B meson 
wave function shape parameter ub = 0.40±0.04 GeV and decay constant Jb = 0.19±0.02 
GeV. The second and dominant theoretical error in these entries arises from the combi- 
nation of the uncertainties of Gegenbauer moments a^aiWh )> a ia ( a 2b )> an d decay 



constant f ai (fbi) m the distribution amplitudes of axial-vector meson Oi (bi). The third 
error is also the combined uncertainty in the CKM matrix elements: p = 0.135looi6 an d 
fj = 0.349to;o}7 [36]. The numerical results implicated that the input parameters(e.g. 
the Gegenbauer coefficients in the nonperturbative wave functions of involved mesons) 
adopted in the pQCD approach should be further improved through the better constraints 
from the experiments to enhance the theoretical precision. 

By comparison, one can easily find that the results on the CP-averaged BRs of B — y 
aidi decays in the pQCD approach agree well with those obtained in the QCDF within 
large theoretical errors, 



BR(B° 

BR(B+ -yofoS) = 21.4; 



— > a{a l 



BR(B° -> a ^) 



54.7+1;^ x 1(T 6 

-12.6 
-10.5 

2.2ii :i 



x 10" 
x 10" 



(In pQCD) 



(71) 



and 



BR(B° -)■ 
BR(B- 
BR(B° -)■ 



a{ a 1 , 
— > afai) 



37.4±i|;| x 10- 6 



22.4+^ 2 3 6 

^ +9 - 3 

u.d_ 0i2 



x 10~ 6 
x 10" 6 



(In QCDF) 



(72) 



in which various errors have been added in quadrature. 

Based on those numerical results given in Tables Ul lII!} some remarks on the CP- 
averaged BRs for B —y a\a\ decays are in order: 



1. As mentioned in the introduction, the experimental measurement for B° —y a; <i 
mode has been performed by BaBar Collaboration and the decay rate is [5|, [lOj 



i "i 



BR(B° -> ata- 



1 a l J Exp. 



47.3; 



-12.2 
-12.2 



X 10' 



(73) 



Combined with the numerical results evaluated with pQCD and QCDF approachesj, 
i.e., Eqs. f JTTj) and f!72|) . one can find that the predicted BR(£?° — y 
well with the present data within uncertainties. 



2 Here, we do not quote BK(B° — » a^a x ) provided in naive factorization [2( for phenomenological 
analysis just because the numerical result is too small to be comparable with the data. 
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TABLE II. Same as Table [J but for B + -»• a+a^bfb^ decays. 



Decay Channels 


B + -+ a+a'j 1 


B+ -> b+b° 


Parameter 


Definition 


This work 


QCDF 


Experiment 


This work 


QCDF 


Experiment 


















BR(10~ 6 ) 


r/r total 


91 , + 7.4+10.0 + 2.2 
z± - -6.5-8.2-1.0 


„„ . + 10.7+6.6 
zz " -8.2-1.5 


- 


- „ + 2.3 + 6.8 + 0.8 
'■ D -1.9-4.1-0.6 


. .+2.5 + 2.8 
*■ — 1. 0-0.0 


- 


fL 


\A L \ 2 


n Q1 +0.00+0.02 + 0.00 
U - M1 ~0. 00-0. 03-0. 00 
n n c+0. 00+0. 02 + 0. 00 
u ' u -0.00-0.01-0.00 

n n q+o. oi+o. 02+0. oo 

u - uo -0. 00-0. 00-0. 00 


74+ 024 
u -'*-0.32 


- 


n c.,1+0.01 + 0.05 + 0.00 
u,o -0. 01-0. 08 -0.00 
n 1n +0. 00+0. 04+0. 00 
u ' iu -0. 01-0. 03-0. 00 
n n7 +0. 00+0. 03+0. 00 
u,u/ -0.01-0.02-0.00 


n nr, + 0.00 
u - ao -0.82 


- 


0||(rad) 
0j_(rad) 

A0i 1 (rad) 
A0 x (rad) 


T + arg ~^- 
f + ar S 33^ 

2 

0j_ — <£x — x 
2 


, 9K +0.02 + 0.08 + 0.00 
^■ za - 0.02 -0.09 -0.00 
o 9 o+0. 02+0. 09+0. 00 
° ° -0.02 -0.10-0. 00 

SB 0.0 

S3 -1.58 


- 


- 


, 7 , + 0.06 + 0-06+0.00 

*-'°-0. 05-0. 05-0. 00 

„ 77 +0. 06+0. 07+0. 00 

-0.05-0.05-0.00 

IK 0.0 

« -1.57 


- 


- 


















.dir 
-A cp 

.dir.L 
-*CP 

,<Hr.ll 

^dir.i 
■^CP 


r-r 
_r+r 
fL —fL 

f_L+fL 

J\\- f u 

/ M +/| 
f±-f± 

f i +/ i 


IB 0.0 
IB 0.0 
is 0.0 
S3 0.0 


" 


- 


ik 0.0 
a: 0.0 
IK 0.0 
IK 0.0 


" 


- 





2. Based on the discussions [3|, |8|, |9[ about the hadron dynamics of axial-vector a\ and 
vector p mesons, due to the similar QCD behavior between them, the decay pattern 
of B — y aiOi modes is therefore analogous to that of B — y pp ones as expected. 



3. Because the decay constant f ai ~ 1.2 x f p w 1.5 x /J, the pattern of BR(i? 



->■ 



ingly. But 
BR(5 



> BR(£> — y pp) is also as expected in the pQCD approach correspond- 



it is worth mentioning that for the color-suppressed decay B — y <x[a 



— > a-iO-iJpQCD 



4 x BR(fi° ->■ 



„o„o\ 

a l a ljQCDF- 



0^0 
1) 



4. Actually, in view of the larger decay constant / Ql and similar QCD behavior of ai 
to p, as a naive estimate, the relation of the BRs among three channels B° — y p°p°, 
- - and B° -> a? a? may be BR(S° -► p°p°) < BR(fi° -> a?p°) < BR(5° -> 



5° -r a^p u 



a^ai). In other words, BR(5 C 



— >■ aVa 



o ^cn 



l "1/ 



-> 



)qcdf should be comparable with BR(5 

BR(fi u -> a?p°) « 2.5 x BR( J B° -> a°a°) in Ref. |3j within the framework of QCDF. 
It is worth stressing that among the above comparisons, the central values of the 
theoretical predictions are adopted for clarification. 



OpQCD- It is thus a bit strange that BR(B° -> p°p°) « 2 x BR(5° -> a^a^) and 



5. Similar to B — y irir and B — > pp decays, one can find that BR(i?° —y afa 1 )> 
BR( J B+ -> atgS) > BR(B° -> a?a$) in the pQCD approach. Moreover, like B° ->■ 
af 7r T decay [37|, in principle, S° — >• a^a]" can be utilized to provide an independent 



measurement of the CKM angle a [6| . Certainly, the currently available statistics is 
too low to perform such an analysis on the angle of a from the current generation 
of B factories. But, potentially, a new generation of Super flavor factories ( Super B) 
are expected to achieve the result with a high luminosity > i0 36 cm~ 2 s -1 |38l. I39J . 

Now let's turn to analyze the phenomenologies of B — y b\bi decays. One can find that 
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the pQCD predictions for BRs of B — y b\b\ decays are 

BR(B° ->■ b+b^) 
BR(5+ -> btb() 
BR(5° -)■ 6?6°) 

and that in the QCDF approach, 



21-4^-1 x 10 -6 ) 

7.61U x 10- 6 , 
29.0tH:\ x 10~ 6 ; 



(In pQCD) 



(74) 



BR(5° -)• 6+6^) 
BR(B+ -> 6+6°) 
BR^ -> b\b{) 



1.0 



+15.8 
-0.8 



X 10" 



1.4±?;[j x 10 
3.2+^ 2 4 3 x 10" 



-6 



(In QCDF) 



(75) 



in which various errors have been added in quadrature too. Though the central values 
of the CP-averaged BRs for B — y b\bi decays in the pQCD approach are much larger 
than those in QCD factorization, they are roughly consistent with each other within large 
errors and in the order of 10~ 6 ~ 10~ 5 . It will be of great interest to measure these 
theoretical predictions to test pQCD and QCDF approaches experimentally. 

From the theoretical predictions presented in Tables Hl lIH[ some discussions on the 
CP-averaged BRs for B — y b\b\ decays are given as follows: 

1. Because of charge conjugation invariance, the decay constant of the l P\ neutral 
meson 6° must be zero. In the isospin limit, the decay constant of the charged b\ 
vanishes due to the fact that the b\ has even G-parity and that the relevant weak 
axial-vector current is odd under G transformation. Hence, f b ± is very small in 

^(0.0028 ± 0.0026)/ fel Q. 



reality. In the numerical calculations, we adopted f b ± 



2. In principle, because of either very small or vanishing decay constant of meson b\, 
the BRs of B — y b\b\ decays should be significantly suppressed by comparing with 
the B — y Oiai ones as naive anticipation. However, as shown in Eqs. (IT5|) - (j22|) . the 
bi( l Pi) meson has the rather different hadron dynamics from its 3 Pi partner, the 



TABLE III. Same as Table UJ but for B c 



a\ a\ , b\ &5 decays. 



Decay Channels 


R<l _. -0 „0 


B° -t- b'l b° 


Parameter 


Definition 


This work 


QCDF 


Experiment 


This work 


QCDF 


Experiment 


















BR(10 -6 ) 


r/r totaJ 


„ „ + 0. 6+1. 6 + 0.1 
z - -0. 5—1. 0-0.1 


„ +0.8 + 9.3 
u -°— 0, 2-0.0 


- 


9Q n+ 7 - 7 + 27 - 9 + 3 - 1 

^ y - -6.9-16.6-2.2 


, 9 + 5.6 + II.O 
°-* — 2.3— 0.8 


- 


Sl 

t\\ 
f± 


\A L \ 2 
\A ± \ 2 


n .,, + 0.01 + 0.02 + 0.00 
u -0.01-0.02-0.01 
n aq+°- 00+0. 00+0. 00 
u ' 4y -0.01-0. 01-0. 00 

n i n+ - 00 + ' 01 + 000 

u,4u -o. 01-0. 02-0. 01 


60+ 000 
u - DU -0-70 


- 


a 1.00 
a 0.00 
S3 0.00 


95+ 002 

u ' yo -0.80 


- 


0||(rad) 
0+(rad) 
A</> 1 (rad) 
A<£^(rad) 


* + arg AZ 
tt + arg -gL. 


„ ro + 0. 05 + 0. 01 + 0. 06 
°-°°-0. 06-0. 92 -0.05 
,,.+0.07+0.02 + 0.06 
J ' D *-0. 07-0. 99-0. 04 

n 77 +0. 10+0.02 + 0.07 
u ' -0.07-1.52-0.09 
n oo + 0. 08 + 0. 02 + 0. 06 
U ' OJ -0. 06-1. 50-0. 08 


" 


- 


. o„+0. 03 + 0. 09 + 0. 04 

^'^ -0.02-0.19-0.00 

. „ .+0.03+0.08 + 0.04 

-0.01-0.15-0.00 

„ „, +0.08 + 0. 16 + 0.09 
u -^_ 0.06- 0.24- 0.00 
. qc + 0. 06 + 0. 13 + 0.07 
*-■ "° -0.05 -0.19-0. 00 


" 


- 


2 

<S> ±_ — 0j_ — ™ 
2 


- 


















^CP 
.dir.L 
■^-CP 

_4dir.i1 

^dir.-L 


r-r 

_r+r 

fL—fL 

lL+fL 

n\~ f \\ 


, 7 „ + 0. 07+0. 00 + 0. 03 
u - ' -0. 04-0. 00-0. 01 

n 90+009+0.05+0.00 

Ul ^ -0.55-0.05-0.01 

„ „ +0.09 + 0.00 + 0.04 
u,yl -0. 04-0. 00 -0.01 

n on+o. 10+0.00+0.04 

U-yu -0. 04-0. 00-0. 01 


" 


- 


n +0.02 + 0.00 + 0.02 
U.Ud_ 01—0 .00-0 .01 
n „, + 0.02 + 0. 00 + 0. 02 
u - u — 0.02— 0.00— 0.01 

„ , Q +0. 00+0. 30 + 0. 00 
u,ly -0. 10-0. 34-0. 15 

n 9, +0.01 + 0. 30 + 0. 00 
u,z -0. 10-0. 33-0. 15 


" 


- 


/ii+/ n 
f±-f± 
/ 1 +/ 


- 
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TABLE IV. The decay amplitudes (in unit of 10~ 3 GeV 3 ) of the B° -)• afa x , bfb x channels 
with three polarizations, where only the central values are quoted for clarification. 



Channel 


B° -)■ af a, 


Decay Amplitudes 


^U 


**L 


A nf, 


A p 


Al fa 


Kfa 


A T 


A p 


L 

N 
T 


3.23 + i8.36 
0.76 + il.98 
1.46 + i3.77 


-0.82 + sO.34 
-0.19 + sO.08 
-0.36 + i0.15 


-0.16 - iO.33 
0.28 + iO.33 
0.57 + 40.71 


0.03 - iO.Ol 
-0.03 + i0.03 
-0.06 + iO.05 


0.47 + i0.09 

0.02 + iO.Ol 

~ 0.00 


-0.11 + bO.29 

-0.00 + sO.01 

~ 0.00 


~ 0.00 

~ 0.00 

-0.02 - iO.OO 


0.31 + i0.21 
-0.41 - s0.61 
0.81 - il.21 


Channel 


B° -> bfb~ 


Decay Amplitudes 


^ 


*f. 


AZn 


A- f . 


< f « 


< 1a 


*■?. 


■*?« 


L 

N 
T 


-0.01 - i0.02 
0.43 + »1.12 
0.81 + i2.09 


~ 0.00 
-0.11 + s0.05 
-0.21 + sO.09 


-1.10 + s3.50 
0.14 - «0.12 
0.32 - i0.21 


-0.36 - sO.09 

0.01 + i0.02 
0.02 + i0.04 


0.52 + il.89 

-0.02 + iO.Ol 

~ 0.00 


-1.26 + s0.42 
~ 0.00 
~ 0.00 


-0.01 - iO.OO 

~ 0.00 
0.01 - iO.OO 


-1.03 - sO.73 
0.17 + i0.13 
0.35 + i0.29 



4. 



ai meson. Correspondingly, one can find the induced anomaly in the theoretical 
pQCD predictions on the CP-averaged BRs and other physical observables shown 
in Tables HHEO 

The B — > &i&i decays receive dominantly large contributions arising from the non- 
factorizable spectator diagrams, which result in the large CP-averaged BRs. To 
clarify this point more clearly, we present the decay amplitudes numerically for ev- 
ery topology with three polarizations in the Tables II VII VI t where only the central 
values are quoted. 

Moreover, the BRs of the B — > bib\ decays exhibit an interesting pattern highly 
different from that of B — > a±ai ones. In terms of the central values in the pQCD 
approach as listed in Tables Ul uTH one can observe that BR(5° — » fe^i) > BR(i?° — >• 



btK) > BR(B+ -> 6+6°), while BR(S° -> 



O^Ch 



ala 



l"l; 



< BR(5- 



— ¥ aja\ 



< BR(B C 



a{ a 1 



Meanwhile, one can also find that BR(B° -»■ &°&°) > BR(B + -> 6+6°) 



— >■ 
> 



BR(B° -»■ bfb^), while BR(B° -> 



o„o> 



a^a 



i w i, 



< BR(B + -> 



a{ a?) 



< 



BR(5° 



— > aj~ % 



in the QCDF approach [3(]. The confirmation of these interesting relations through 
the relevant experiments may shed light on the QCD dynamics involved in these 
considered channels. 



Although the components of B — V b\b\ decays at the quark level are same as that of 
B — > pp and B — » a\a\ decays, the different phenomenologies have been shown(See 
Tables Hl lHip because of the different QCD behavior between cq and b\ mesons. One 
can therefore expect that B — > b\b\ decays will imply some new information on the 
CKM unitary angle, reliability of pQCD approach, and so on. 



As mentioned in the above, our pQCD prediction on the BR of B° — > o + o5~ decay 
is consistent with the data reported by BaBar Collaboration very recently. Though the 
theoretical errors are a bit large, the central value of this channel can still be employed 
to roughly estimate the BRs of other decay modes. Here we define four parameters, say, 



» R 0H 



K bb : 



»-r00» 
K afe 



and "R 



ab 



as follows, 



R„ 



t b+ BR(B° -> a 



i "i 



r B o BR(B^ 



->■ a{a\) 



2.7, 



R°+ = 



t b+ BR(B° -j- 6+67/) 

T^o" ' ERjm~Tb[bJ) 



3.0 ; (76) 
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TABLE V. The decay amplitudes (in unit of 1CT 3 GeV 3 ) of the B + -)■ a\a^ 
three polarizations, where only the central values are quoted for clarification. 



bfb® channels with 



Channel 


S+ -> a+a? 


Decay Amplitudes 


^s 


*?. 


Aln 




*z f . 


<f* 


* T t. 


*% 


L 

N 
T 


1.98 + »5.13 
0.48 + il.25 
0.92 + i2.37 


-0.18 + i0.07 
-0.04 + J0.02 
-0.08 + W.03 


0.18 + i0.34 
-0.25 - sO.27 
-0.49 - s0.60 


-0.01 + sO.OO 
0.01 - iO.Ol 
0.01 - iO.Ol 


- 


- 


- 


- 


Channel 


B+ -H.+ 6? 


Decay Amplitudes 


*?. 


*?. 


< fa 


<r. 


<f» 


<f» 


A % 


*% 


L 

N 
T 


-0.01 - 40.01 
0.26 + i0.68 
0.49 + il.27 


~ 0.00 
-0.02 + iO.Ol 
-0.04 + i0.02 


0.59 - J3.30 
-0.09 + i0.10 
-0.20 + i0.20 


0.11 - i0.02 

~ 0.00 
-0.01 - J0.01 


- 


- 


- 


- 



BR(B° -»■ a+o^) 
"'' " BR(B° ->■ b+b^) 



-ruu _ 



2.6 



R ab + = 



BR{B + -»■ ata\) 
BR(5+ ->■ 6+6?) 



2.8 



(77) 



We expect the above four ratios could be tested at the ongoing LHC and forthcoming 
Super-B experiments. 

Finally, we should stress that both of the color-suppressed modes B° — > a®a® and 



1"! 



B° —7- b\b\ themselves exhibit the dramatically different phenomenologies. Similar to 
B° —7- 7r°7r°, Tc°p ,p°p decayed, B° — > ct^a® and b®b® decays are closely related to the 
color-suppressed tree amplitudes C u oc a 2 (= C\ + C2/3), where C\ and C 2 are Wilson 
coefficients. At leading order, the sign of C 2 is positive while the sign of C\ is negative, 
which can cancel each other mostly. For example, the numerical result of a 2 is about 
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Furthermore, one can 
and AT f „ from factorizable 



1.1 x 10~ 3 when the running hard scale is taken at \x = 2.5 GeV 

easily find from the Table [VI] that the tree contributions A^ s a±±u. s\ n t s 

spectator and nonfactorizable spectator diagrams cancel each other significantly for B° — > 

a\a\ channel. Thus BR(5° — > a^a^) is rather small relative to BR(5° — > a^a^) and 

BR(B + — > a + o°) q It should be stressed that the small quantity a 2 in B° — > b^b® need 

not to be considered seriously because the decay constant f b o is exact zero. But, the 

resulting BR(5° — > £>5^i ) * s SUCI1 large that reaching 29.0 x 10~ 6 numerically. It will be 

highly interesting to measure this rate to test the availability of pQCD approach in the 

channels with X P\ mesons. 



3 As for the color-suppressed processes in the decays of B mesons, which have been extensively studied in 
plenties of literatures with various of methods and/or schemes within and beyond the standard model. 
However, to our best knowledge, they seem to be a longstanding "puzzle" in B physics because one 
can not resolve it self-consistently in the current approaches/methods. 

Unfortunately, up to now, the color-suppressed tree amplitude C seems to be an important but the 
least understood quantity in B meson decays [24 1 

5 Recently, the authors in Ref. [24j proposed a solution to the B —> nn puzzle by considering the 
contributions arising from " Glauber-gluon- region" . However, it is worth mentioning that this class of 
contributions may make very little effects to the results on the considered B — > a\ a?, b\ b\ decays in 
the present work, which because, as argued in the literature, the Glauber effects can only contribute 
significantly to the pion but much less to the p meson. 
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TABLE VI. The decay amplitudes (in unit of 10" 3 GeV 3 ) of the B c 



a? a 



i , b\ b\ channels 



with three polarizations, where only the central values are quoted for clarification. 



Channel 


H° _1 r.° „« 


Decay Amplitudes 


^ 


^ 


Al u 


<U 


K 1a 


< fa 


*?. 


*?. 


L 

N 
T 


0.30 + 40.79 
0.06 + 40.15 
0.12 + 40.30 


-0.40 + 40.17 
-0.09 + 40.04 
-0.17 + 40.07 


-0.28 - 40.58 
0.45 + 40.51 
0.88 + 41.10 


0.03 - 40.01 
-0.03 + 40.02 
-0.06 + 40.05 


0.31 + 40.07 

0.02 + 40.01 

~ 0.00 


-0.08 + 40.20 
~ 0.00 
~ 0.00 


~ 0.00 

~ 0.00 

-0.01 - 40.00 


0.22 + 40.14 
-0.29 - 40.43 
-0.59 - 40.86 


Channel 


B° -+ b'( b'l 


Decay Amplitudes 


**t. 


^ 


Ku 


<U 


K fa 


<f. 


^ 


*?* 


L 

N 
T 


0.00 
0.04 + 40.11 
0.08 + 40.21 


0.00 
-0.05 + 40.02 
-0.10 + 40.04 


-1.34 + 45.77 
0.19 - 40.18 
0.42 - 40.35 


-0.30 - 40.12 
0.01 + 40.01 
0.02 + 40.03 


0.34 + 41.32 

-0.01 + 40.01 

~ 0.00 


-0.90 + 40.31 

~ 0.00 
~ 0.00 


~ 0.00 

~ 0.00 

0.01 - 40.00 


-0.73 - 40.52 
0.12 + 40.10 
0.25 + 40.20 



B. Polarization Fractions 



We have also computed the polarization fractions for B — > a^i and &1&1 decay modes. 
Based on the helicity amplitudes ( 12^1) . we can define the transversity amplitudes as 



Al = -Zm%M L , A\ = tV2m 2 B M N , A ± = ^r 2 r 3 y/2(r 2 - l)m 2 B M T ■ (78) 

for the longitudinal, parallel, and perpendicular polarizations, respectively, with the nor- 
malization factor £ = ^jG 2 F P c /{lQ r Km B Y) and the ratio r = P 2 • P 3 /{m 2 B r 2 r3). These 
amplitudes satisfy the relation, 



|.4 L | 2 + |-4|j| 2 + |-4j 



1 



(79) 



following the summation in Eq. flTOj) . The polarization fractions fi,f\\ and f± can thus 
be read as, 



fi 



,-L) 



\A. 



H\U)\ 



\A L \ 2 + \A\\\ 2 + \Ai 



\A 



L(\U)\ > 



?0) 



The numerical results of fractions with three polarizations for B — >■ a\a\ and B — > b\b\ 
decays in the pQCD approach have been presented in Tables HllIIII Based on these values, 
we give some phenomenological analysis: 



1. As discussed in Sec. IIVAI the prediction on BR(_B° — > afa~[) with pQCD and QCDF 
approach, respectively, is in good agreement with the measurement given by BaBar 
collaboration. However, the fraction of longitudinal polarization for B° — > afai is 
not the case. Theoretically, this considered channel is dominated by the longitudinal 
contributions, 



/ L ( J B°^a+ar) pQC D = 0.76] 
/ L ( J B°^a+ar) QC DF = 0.64] 



0.03 

-0.04 > 

-0.07 . 

■0.17 ) 



51) 

J2) 



and fi predicted in pQCD and QCDF are very close to each other; on the other 
hand, experimentally, 



f L (B° -)> af oJexp. 



U - O1 -0.24 



J3) 
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it seems to be governed by the transverse ones. But, it should be mentioned that 
the measurement performed by BaBar collaboration still have very large errors and 
should be greatly improved, in order to test the theoretical predictions in the near 
future. 

2. Because the QCD behavior of axial-vector a\ meson and that of vector p meson are 
analogous to each other, the numerical results show the similar pattern of longitu- 
dinal polarization between B — >■ oi Oi and B — Y pp 33|, |34J decays: B° — > a 1 af and 



B ± — > a x a® decays are dominated by the longitudinal component, reaching around 
76% and 91%, respectively, while B° — > a® a® is governed by the transverse one at 
leading order, fi « 12%. Maybe the situation for B° — > a\a\ will be highly im- 
proved after taking the higher QCD corrections into account. (See B° — > p°p° [32| 
for example.) 

3. Similar to B° — > a^a~[ and B + — > afa^ modes, B° — > b~l 6f and B + — > b\b\ ones are 
also dominated by the longitudinal polarization components. However, dramatically 
different from B° — > a\ a!{ decay, B° — > b\ b\ channel is absolutely governed by the 
longitudinal polarization contributions, say, fi ~ 100%. 

4. For a clear clarification to the polarization fractions, we have listed the contributions 
from every topology with three polarizations in the Tables IIVHVII Generally speak- 
ing, from Tables HVtiVIt one can find that relative to B — > a\a\ decays, all three 
B —7- b\b\ decays suffer from large nonfactorizable spectator tree contributions in 
the longitudinal component, which therefore lead to the dominance of longitudinal 
polarization fraction. In contrast to B° —¥ 6^1 decay, B° —$■ a^ai channel receives 
a bit large contributions from nonfactorizable spectator tree diagrams in both of 
transverse polarizations. 

We expect the above observations would be tested by the future experiments, then 
could provide more information for understanding the underlying helicity structure in 
these types of decays. 



C. Effects of Nonfactorizable Spectator and Annihilation Contributions 

To see whether weak annihilation contributions play important role in these considered 
decay modes, we test the CP-averaged BRs and longitudinal polarization fraction by 
neglecting the annihilation diagrams, which can not be perturbatively calculated in the 
QCDF approach. Moreover, as claimed in the references within the framework of QCDF, 
the hard spectator scattering contributions also suffer from endpoint singularities at the 
level of twist-3. In other words, the calculations with these terms in the QCDF need 
always the adjustments based on the measurements at relevant experiments. 

Without the annihilation contributions in both B° — > a^ai and B° — > a^cii decays, we 
find the following CP-averaged BRs and longitudinal polarization fractions through the 
numerical evaluations in the pQCD approach, 

BRt^B -> a^O w 51.0 x 10~ 6 , f L (B° ->■ a+a^) » 0.77 ; (84) 

BR(B° ->■ a ^) w 1.5 x 10" 6 , f L (B° -> a ^) w 0.07 . (85) 
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which means that the annihilation contributions account for a small ratio and could be 
neglected safely in B° — > afa^ mode, while it is not the case in the decay of B° — > a^a®. 
When we neglect the decay amplitudes arising from both of nonfactorizable spectator 
diagrams and annihilation ones, the numerical results for the CP-averaged BRs and lon- 
gitudinal polarization fractions of B — > a\a\ decays are as follows: 

BR(B° -* a+a^) « 49.6 x 10~ 6 , f L (B° -► a+a^) « 0.84 ; (86) 

BR(B + -► a+a\) « 20.1 x 10~ 6 , f L (B° -)- a+a\) « 0.84 ; (87) 

BR(£° -> o?a$) « 0.5 x 10" 6 , f L (B° -> a?a$) « 0.88 . (88) 

One can find that the nonfactorizable spectator contributions in the B° — > afa± and 
£> + — > afa® modes are so small that they could be neglected, while that in B° — >■ a^a? is 
large and play an important role. This phenomenon can be easily found from the decay 
amplitudes shown in Tables HVllVIl In the above Eqs. (184]) - ( 188]) . only the central values 
are quoted for clarification. 

Similarly, we predict the CP-averaged BRs and longitudinal polarization fractions with- 
out the factorizable and nonfactorizable annihilation diagrams for B° — > bfbi, b\b\ decays, 

BR(B° -)> 6+6^) « 9.5 x 10~ 6 , f L (B° -)■ bfb^) « 0.74 ; (89) 

BR(5° -> 6°6°) « 18.7 x 10" 6 , / L (£° -> 6°6?) « 0.99 . (90) 

which means that there should exist large contributions from annihilation diagrams in 
these two decays. Meanwhile, they also indicate that there are large nonfactorizable 
spectator diagrams [l[ due to the fact of large longitudinal polarization fractions and 
extremely small or vanished decay constant in the longitudinal twist-2 wave function. 
By considering only factorizable emission diagrams in Fig. [1] the predicted BRs in the 
pQCD approach are determined completely by the transverse components because the 
longitudinal contributions from Fig. 1(a) and 1(b) are sharply suppressed by the tiny or 
zero decay constant, 



BR(B° -► 6+6^) « 2.6 x 10~ 6 , 


f L (B°^btK)z 


rfO.O; 


(91) 


BR(B + ->■ 6+6°) « 1.0 x 10~ 6 , 


f L (B° -». 6+6?) p 


^0.0; 


(92) 


BR(B° ->■ 6°6°) « 3.0 x 10~ 8 , 


/ L (5° -► 6?6?) c 


a 0.0. 


(93) 



which exhibit evidently the dominated spectator and/or annihilation contributions in- 
volved in the B — > b\bi channels. Notice that for B° -> 6565 the CP-averaged branching 
ratio is significantly small just because the emission decay amplitudes are proportional to 
the combined Wilson coefficient a2 in both transverse polarizations. 

D. Direct CP- violating Asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries for B — > a%ai and 
b\b\ decays in the pQCD approach. It is conventional to combine the three polariza- 
tion fractions in Eq. (180]) with those of its CP-conjugate B decay, and to quote the six 
resulting observables corresponding to tranversity amplitudes as direct induced CP asym- 
metrie^f • 



(i 



The direct CP asymmetries in transversity basis can be defined as 

^ = 7^> = £,||,1_) (94) 

where the definition of / is same as that in Eg. ([80]) but for the corresponding B decay. 
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As for the direct CP-violating asymmetry in these considered modes, considering the 
involved three polarizations, whose definitions are as follows, 

,dir s r-r = |M(^->/ final )| 2 -|M(i?^/ final )| 2 
CP r + r \M{B^J^)\* + \M{B^h^)\^ l } 

where Y and M denote the decay rate and decay amplitude of B — > a\a\,b\b\ decays, 
respectively, and T and M are the charge conjugation one correspondingly. 

Based on the above definitions on direct CP-violating asymmetry and numerical cal- 
culations in the pQCD approach(see Tables Hl lIIII) . some remarks are as follows: 

1. The direct CP asymmetries for B — > c^ai decays in the pQCD approach can be 
read as, 

Acp(B° ->• afaf) « -0.04 ± 0.01 , (96) 

A d j P {B ± -* afa\) « 0.00 , (97) 

Af P (B° -> a?o?) « -0.78+°;°* . (98) 



which are very similar to those in the B — > pp decays [34J correspondingly, where 
the various errors as specified have been added in quadrature. 

2. As for the direct CP-violating asymmetries of B — > b\ b\ decays, 

Af P (B° -► bfbf) « -0.00+H1 , (99) 

^(5 ± ->■ 6f 6?) « 0.00 , (100) 

^(S° ->• 6?6?) « -0.03+^1 • (101) 

One can find that the numerical results in the pQCD approach at leading order are 
very small, even to be zero within uncertainties as presented in Tables Hl lIIII 

3. Meanwhile, we calculate the direct CP-violating asymmetries in every polarization 
and give the results in the pQCD approach as 

/jdir.Z, _ n ,,+0.02 Adit,\\ _ _ n pro+0.13 4dir,± _ _q ^+0.14 . Mn.cy\ 

•^CP ~ u - ±± -0.02 i ^-CP ~ u - oz -0.10 ' ^CP ~ u - d ^-0.09 i l iuz J 

for B° — )■ a 1 af mode, and 

j.dir,L _ n on+o-io ^ dir 'H - _nqi+o.io ^dir,_L _ _ n on+o.n . / lf) o\ 

J\ CP — U.ZU_ 055 , ^-CP ~ u - J± -0.04 5 ~^-CP ~ u - yu -0.04 i \ LVO ) 

for B° — >• a^ai channel, and 

yidir.L _ _ n nf -+0.04 Adir,\\ _ n q/|+0.05 /|dir,± _ n oo+0.07 . flfl/H 

•^CP ~ u - ud -0.05 > ^CP ~ u -°^-0.08 J ^CP ~ u -°°-0.07 ' \ LV> ^> 

for B° — > b 1 bf mode, and 

A dir,L _ _ n n o+0.03 yidir,|| _ n i q+0.30 jdir,± _ n 9 o+0.30 . / ln c\ 

J\ CP — U.UO_ 00 2 , -^-CP ~ u - ±y -0.38 J -^CP ~ U - ZO -0.38 i \ lvo ) 

for B° — > b\b\ channel, respectively, in which the various errors as specified have 
also been added in quadrature. These direct CP-violating asymmetries are expected 
to be confronted with the relevant measurements in the future. 
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4. Because of the lack of strong phase from the annihilation diagrams and the rather 
negligible contributions just from electroweak penguin operators in nonfactorizable 
spectator diagrams, the direct CP violations in the B ± — y afa® and bfb® decays is 
absent for every polarization naturally, which can be seen easily in Table |V] 

We also define another two quantities to reflect the existence of direct CP-violating 
asymmetries indirectly, 

A0|| = and A0 ± = , (106) 

where </>n and 0j_ are the CP-conjugated ones of relative phases <fi\\ and <p±, respectively. 
Based on the definitions of transversity amplitudes, the relative phases 4>\\ and <j>±_ are 
defined as, 

0|| = arg — i and (j)± = arg — - , (107) 

The theoretical predictions of relative phases for B — y a\ a\ and &i b\ modes in the pQCD 
approach have been presented in Tables Hl lIIIl which will be tested by the measurements 
at B factories, ongoing LHC, even forthcoming Super-B experiments. Note that the 
definitions of Al,||,_l as given in Eq. (1751) are consistent with those in [3], except for an 
additional minus sign in Al, so that our definitions of the relative strong phases <j>\\,±. (see 
Tables Hl lHip also differ from the ones in [3| by it, which is added to cancel the additional 
minus sign in the definition of Al in Eq. 



At last, it is worth of stressing that the theoretical predictions in the pQCD approach 
still have large theoretical errors(See Table [J for example) mainly induced by the still 
large uncertainties of distribution amplitudes from the shape parameter ujb of heavy B 
meson and the Gegenbauer moments a,- of light axial- vector oi and b\ mesons. We 
need the nonperturbative QCD efforts and experimental constraints to effectively reduce 
the errors of these essential inputs. Any progress at this aspect will help us to improve 
the precision of the pQCD predictions. 



V. SUMMARY 

In this work, we studied the charmless hadronic B — y a\a\ and b\b\ decays by employ- 
ing the pQCD approach based on the kr factorization theorem. We calculated not only 
the factorizable emission diagrams, but also the nonfactorizable spectator and annihila- 
tion ones. Our theoretical predictions in the pQCD approach will provide an important 
platform for testing the SM and exploring the helicity structure of these considered decays 
and the hadronic dynamics of the axial-vector a\ and b\ mesons. They can also provide 
more information on measuring the unitary CKM angles and understanding the decay 
mechanism of color-suppressed modes. 

The pQCD predictions for B —y OiOi, b%bi channels are displayed in Tables (Ifl lHip . 
From our evaluations and phenomenological analysis, we found the following results: 
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The CP-averaged branching ratio of B° — > afa^ mode in the pQCD approach is in 
good consistency with that given by preliminary measurement and that presented 
in the QCDF framework, respectively, within errors. 

The pQCD predictions for the CP-averaged branching ratios of B — > a%a\, b\b\ 
decays are in the range of 10~ 5 to 10~ 6 , which can be easily accessed at the B 
factories of BaBar and Belle, running LHC, and forthcoming Super-B experiments. 

The numerical results in the pQCD approach, specifically, on the CP-averaged 
branching ratios and longitudinal polarization fractions of the considered B — > 
ciidi, bibi decays are basically consistent with those given in the QCDF framework, 
except for /l(-B° — > 



cucn 



a\a 



• The theoretical predictions in the pQCD approach have large uncertainties, which 
mainly arise from the nonperturbative input parameters with still large errors, for 
example, the distribution amplitudes describing the hadron dynamics of the involved 
mesons. We expect these inputs will be well constrained when more data become 
available. 

• We here simply take the short-distance contributions into account in the evaluations 
of B — > OiOi, 61&1 decays. Maybe the final state interactions for these considered 
modes play an important role, more relevant studies are therefore helpful for us to 
provide reliable pQCD predictions. 
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